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SPACE-FILLING CURVES OF SELF-SIMILAR SETS (I): ITERATED 
FUNCTION SYSTEMS WITH ORDER STRUCTURES 


HUI RAO AND SHU-QIN ZHANGf 

Abstract. This paper is the hrst paper of three papers in a series, which intend to 
provide a systematic treatment for the space-filling curves of self-similar sets. 

In the present paper, we introduce a notion of linear graph-directed IFS (linear GIFS 
in short). We show that to construct a space-filling curve of a self-similar set, it is amount 
to explore its linear GIFS structures. Some other notions, such as chain condition, path- 
on-lattice IFS, and visualizations of space-filling curves are also concerned. 

In sequential papers [7] and [^, we obtain a universal algorithm to construct space- 
Hlling curves of self-similar sets of finite type, that is, as soon as the IFS is given, the 
computer will do everything automatically. Our study extends almost all the known 
results on space-filling curves. 

MSG 2000: 28A80, 37A05,37B10. 


1. Introduction 

Space-filling curves have fascinated mathematicians for over a century. Its history 
started with the monumental result of Peano in 1890 ([ 22 ])- One year later, Hilbert 
gave an alternative construction, now called Hilbert curve. In 1921, Sierpihski discovered 
Sierpihski space-filling curve, and it was generalized by Polya. (See [5].) For variations of 
the above constructions, see the survey book of Sagan [25] . 



Figure 1. Space-filling curves of Peano, Hilbert and Sierpihski. 
t The correspondence author. 
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Around 1970’s, several remarkable progresses have been made: J. Heighway, a physicist, 
found the Heighway dragon ([IS [8]); W. Gosper, a computer scientist, found the Gosper 
island m); H. Lindenmayer, a biologist, introduced L-system (IZl), which becomes a 
powerful method to produce space-filling curves later. 




Figure 2. The Gosper island is a 7-reptile, and the Gosper curve is an 
optimal parametrization. 

Two important facts are gradually recognized: All the constructions are based on cer¬ 
tain self-similar structures, and certain ‘substitution rules’ play an essential role in the 
constructions. 

Next major progress was made by Dekking [9] (1982), where he claimed that he “in¬ 
troduce a powerful method of describing and generating space-filling curves”. This paper 
has important impact on both space-filling curves and fractal geometry. On the fractal 
geometry aspect, leads to the emerge of the notion of graph-directed iterated func¬ 
tion system. On the space-filling curve aspect, Dekking’s method has been accepted by 
computer scientists and as the “vector method”. 

In recent years, various interesting constructions of space-filling curves appear on the 
internet, for example, “www.fractalcurves.com” (see m) and “teachoutl.net/village/” 
(see [27]). Besides, space-filling curves of higher dimensional cubes have been studied by 
Milne [21] and Gilbert [Uj. For applications of space-filling curves, see Bader |3| and the 
references therein. 

In this paper and two sequential papers, we unveil the mystery of space-filling curves 
by providing a rigorous and systematic treatment. 

First, let us specify our meaning of space-filling curves. We call an onto mapping 
from an interval [a, b] to a self-similar set K an optimal parametrization, if it is almost 
one-to-one, measure-preserving and 1/s-Holder continuous, where s = dinij;/AT is the 
Hausdorff dimension of K. (For precise definition, see Section [2l) It is observed that most 
classical space-filling curves fulfill the above requirements mm), while some others 
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Figure 3. Heighway dragon curve and Lebesgue curve. 

like the Lebesgue curve does not (see Figure [3l( right)). It is proper to call an optimal 
parametrization a space-filling curve if K has non-empty interior, and call it a fractal-filling 
curve otherwise. However, for simplicity, we shall just call an optimal parametrization a 
space-filling curve. 

The main contribution of this paper is that we introduce a notion of linear GIFS to 
describe and handle space-filling curves. The graph-directed iterated function system, or 
GIFS in short, is an important notion in fractal geometry. We equip the functions in a 
GIFS with a partial order and call it an ordered GIFS, and this order induces a dictionary 
order of the associated symbolic space. An ordered GIFS is called a linear GIFS, if 
every two consecutive cylinders have non-empty intersections (see Section [3] for precise 
definition). We show that 

Theorem 1.1. Let he the invariant sets of a linear graph-directed IFS satisfying 

the open set condition and 0 < < oo for j = 1,... ,N, where 6 is the similarity 

dimension, then Ej admits optimal parametrizations for every j = 1,... ,N. 

The proof of Theorem 11.11 is constructive; hence, to construct space-Hlling curves is 
amount to seek a linear GIFS structure of the given set. The common point of the L- 
language method and Dekking’s vector method is that, first they construct a linear GIFS, 
and then verify the open set condition. 

Remark 1.1. (i) The notion of linear GIFS can be regarded as a completion of the study 
of Dekking [S] . 

(ii) After we finished this paper, we acknowledge that an idea similar to our linear 
GIFS has appeared in Akiyama and Loridant mm when studying the parameterizations 
of boundaries of self-afRne tiles. 
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For an ordered GIFS, one can associate to each invariant set Ej a head (the point with 
the lowest coding) and a tail (the point with the highest coding). Using heads and tails, 
we define a chain condition (see Section 01) which provides a simple and practical criterion 
of linear GIFS. 

Theorem 1.2. An ordered GIFS is a linear GIFS if and only if it satisfies the chain 
condition. 

To ‘see’ a space-filling curve, we need to visualize or to approximate a space-filling 
curve. Using linear GIFS, in Section [6l we give a precise definition of visualizations of a 
space-filling curve. 

To illustrate our theory, we give a brief introduction to the path-on-lattice IFS in Section 
[5j A nice collection of space-filling curves given by path-on-lattice IFS, many of them are 
well-known, can be found in the website [29]. A detailed study of the path-on-lattice IFS 
can be found in [25] . 

To find the linear GIFS structure of a given self-similar set is a hard question. This 
question is studied in sequential papers [7] and [23] . We show that 

Theorem 1.3. f[7] and m) Let K he a connected self-similar set satisfying the open set 
condition, if K has the finite skeleton property, then it admits optimal parametrizations. 
In particular, if K satisfies a finite type condition (another important condition in fractal 
geometry), then it possesses finite skeletons and hence admits optimal parameterizations. 

Our theory gives a universal algorithm to find space-filling curves of self-similar set of 
finite type, that is, as soon as the IFS is given, the computer will do everything. Our 
study extends almost all the known results on space-filling curves, and shows the internal 
relation between the space-filling curve and the recent developments of fractal geometry. 

Example 1.2. The four-star tile. Pictures in Figure 0] are taking from [27], but there is 
no explanation how to obtain the space-filling curve. Our study will fill all the gaps from 
the left picture to the right in Figured] which is interesting and highly non-trivial ([7]). 
A sketch of the approach is provided in Section |7J 

The paper is organized as follows. In Section [21 we define optimal parametrization for 
general compact sets. We introduce the linear GIFS and the chain condition in Section [3| 
and Section 0] respectively. Section [5] is devoted to the path-on-lattice IFS on the plane. 
Visualizations of space-filling curves are discussed in Section [6l Section [7] studies the 
four-star tile. In Section [S] we prove Theorem ll.il using a measure-recording GIFS. 
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Figure 4. The four-star tile and a space-filling curve. 

2. Optimal parameterizations of self-similar sets 

Let K C be a non-empty compact set. We call K a self-similar set, if it is a union 
of small copies of itself, precisely, there exist similitudes S\,... ,Sn : —)• such that 

N 

K = y Sj(K). 

i=i 

In fractal geometry, the family {S'!,..., ^at} is called an iterated function system, or IFS in 
short; K is called the invariant set of the IFS mm- We denote by the s-dimensional 
Hausdorff measure. A set E C is called an s-set, if 0 < < oo for some s > 0. 

The IFS {S'!,..., Sn} is said to satisfy the open set condition (OSC), if there is an open 
set U such that Ui^i Si{U) C U and the sets Si{U) are disjoint. It is well-known that, if 
a self-similar set K satisfies the open set condition, then it is an s-set. (See [llj.i 

Remark 2.1. If an IFS satishes the OSC condition, and dim// K equals the space dimen¬ 
sion, then K has non-empty interior ([2^), and it is a self-similar tile. Especially, if the 
contraction ratios of Si are all equal to r, then K is called a reptile. (In this case, we must 
have r = 1/v^, where d is the dimension of the space.) 

Motivated by the studies of the space-filling curves, it is natural to define an optimal 
parametrization of more general sets (see m)- Denote £ the one-dimensional Lebesgue 
measure. 

Definition 2.2. Let K C be an s-set. An onto mapping : [0,1] —> AT is called an 
optimal parametrization of K if the following three conditions are fulfilled. 

(z) -ip is almost one-to-one, precisely, there exist K' <Z K and I' C [0,1] such that 
71^{K \ K') = £([0,1] \ /') = 0 and ip : ^ AT' is a bijection; 
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(ii) ijj is measure-preserving in the sense that 


= cC{F) and C{il^-^{B)) = 


for any Borel set F C [0,1] and any Borel set B G K, where c = TL^{K). 
{in) ^|J is 1/s-Holder continuous, that is, there is a constant c' > 0 such that 

|'0(x) — '0(y)l < c'|x — for all x,y £ [0,1]. 


Our main concern is: Does every connected self-similar set admit an optimal parametriza- 


tion? According to the theorem of Mazurkiewicz-Hahn ([25]), a set is the image of [0,1] 
under a continuous mapping if and only if it is compact, connected, and locally connected. 


We note that a connected self-similar set fulfills these conditions, since a self-similar set is 
locally connected as soon as it is connected f |15]l. 



Figure 5. Optimal parameterizations of the Sierpihski gasket and carpet. 

A detailed study of the carpet is carried out in [7] . 

For parameterizations of fractal sets, the previous studies focused on the Holder con¬ 
tinuity. de Rham uni, Hata m and Remes [23] showed the existence of 1/s-Holder 
continuous parameterizations for certain classes of self-similar sets. Akiyama and Lori- 
dant mm showed the existence when K is the boundary of a class of self-affine tiles (their 
motivation is to provide an alternative way to show the disk-like property of some planar 
tiles). Martin and Mattila [19] gave some negative results when K is disconnected. 


3. Linear GIFS 


In this section, we introduce the notion of linear GIFS. 

Let us start with the definition of GIFS. Let G = {A, F) be a directed graph with vertex 


set A and edge set F. Let 
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be a family of similitudes. We call the triple (^, F, Q), or simply G, a graph-directed iterated 
function system (GIFS). We call (^, F) the base graph of the GIFS. Very often but not 
always, we set .A to be {1,... , N}. 

Let Tij be the set of edges from state i to j. It is well known that there exist unique 
non-empty compact sets {Ei}^^ satisfying 

N 

(3.1) = U U 9~({Ej), l<i<N. 

We call {Ej}^^^ the invariant sets of the GIFS ('[20][3]L 

We say the above GIFS satisfies the open set condition (OSC), if there exist open sets 
Ui ,..., Un such that 

N 

U U g^{Uj) cUi, l<i<N, 

j=i -yeTij 

and the left-hand sides are non-overlapping unions ([ 20 ] [ 12 ]). 

Remark 3.1. It is seen that the set equations (13.1|) give all the information of a GIFS, 
and hence provide an alternative way to define a GIFS. We shall call (13.ip the set equation 
form of a GIFS. 

3.1. Symbolic space related to a graph G. Let G be a directed-graph. A sequence 
of edges in G, denoted by cj = a;ia ;2 • ■ ■ ^n, is called a path, if the terminate state of tVi 
coincides with the initial state of cuj+i for 1 < i < n —1. We will use the following notations 
to specify the sets of finite or infinite paths on G = (A., F). For i G A, let 

Ff, F: and F- 

be the set of all paths with length k, the set of all paths with finite length, and the set of 
all infinite paths, emanating from the state i, respectively. Note that Ft = 

For a sequence cj = {tvk)^=i, set cj|„ = cJitJ 2 ... cj„ be the prefix of cj of length n. For 
an infinite path cj = {cOn)^=i G F^, we call 

[cJi . . . CJfi] .— {'y G Fj , 'y\n — • • • ^n} 

the cylinder associated with cji... cj„. 

For a path 7 = 71 ... we denote 


E-^ ■ g-fi o • • • o (F/j(..y)), 
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where t{'y) denotes the terminate state of the path 7 (also 7 „). Iterating (1,3.11) fc-times, 
we obtain 

(3.2) U 

We define a projection vr : (rf°,..., T^) —..., where vr* : is defined 

by 

(3.3) {vr,(a;)} := f| 

n>l 

For X G Ei, we call u a coding of a: if 7 rj(cj) = x. It is folklore that 7 rj(F^) = Ei. 

3.2. Order GIFS and linear GIFS. Let {A,T^Q) be a GIFS. To study the ‘advanced’ 
connectivity property of the invariant sets, we equip a partial order on the edge set F 
enlightened by set equation (13.2p . Let Fj = F)^ be the set of edges emanating from the 
vertex i. 

Definition 3.2. We call the quadruple (.4., F,^,-<) an ordered GIFS, if 7 is a partial 
order on F such that 

(i) 7 is a linear order when restricted on Fj for every j G A] 

{ii) elements in Fj and Fj are not comparable li i A 3- 

We denote the edges in Fj by 7 ^, 1 , 71 , 2 , • • • AiA ™ ascending order. For simplicity, 
we use the following equations to describe an order GIFS (equation form of an ordered 
GIFS): 

— 07^1 (■®i(7i,l)) + ■ ■ ■ + ~ 1) • • • ) 

where we use ‘+’ instead of ‘U’ to emphasize the order. 

The order 7 induces a dictionary order on each F^, namely, 7172 ... 7 ^ a;ia ;2 ... if 
and only if 71 ... 7 £_i = wi... ooi-i and 7 ^ -< for some 1 < I <k. Observe that (F^, -<) 
is a linear order. Now we can define the linear GIFS. 

Definition 3.3. Let {A,T,Q,-<) be an ordered GIFS with invariant sets It is 

termed a linear GIFS, if for all i G .4. and fc > 1, 

E~f n E^^ A 0 

provided 7 and cj are adjacent paths in F^. 
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3.3. Linear IFS. An IFS {S'!,..., S’at} is a special class of GIFS, where the vertex set is 
a singleton, and the edge set consists of N self edges which we denote by 1,..., A^. The 
IFS becomes an ordered IFS, if we assume the natural order 1 ^ 2 -<•••-< A^. 

Remark 3.4. For an IFS, the associated symbolic space is much simpler. Let S = 
{!,..., A^}. For m > 1, we denote = {1,2- S* = and = 

{1, 2 • • • A^}^. For convention, instead of calling uji.. .ujn £ S"' a path, we call it a word. 
For iii 2 .. .im ^ S*, we call {iii 2 ■ ■ ■ im] = {^ S F,°° ; u\m = iii 2 ■ ■ ■ im} a cylinder. 

Denote by K the invariant set of the IFS. The projection map vr : K is 

OO 

1 

where = ‘S'n o Sig o • • • o Si^. We call u a coding of x if 7 r(aj) = x. 

Clearly, the von Koch curve is generated by a linear IFS. In Section O we shall show 
that the Peano curve is generated by a linear IFS, while the Heighway dragon and the 
Hilbert curve are generated by linear GIFS’. 

Example 3.5. Sierpiriski curve. Let Ti,T 2 ,r 3 and be a partition of the unit square 
indicated by Figure [ 6 )^right). In Figure [ 6 )^left), Tj, 1 < j < 4, is divided into 4 small 
triangles, where the numbers indicate the order of the small triangles. The four small 
triangles are images of Tj under a map of the form (x + b)/2. Hence, we obtain a linear 
GIFS. Precisely, 

_ Ti T2 Ti + 1 Ti + 1 

2 2 2 2 

Similarly equations can be obtained for T 2 ,T 3 and r 4 . 



Figure 6 . A linear GIFS which generates the Sierpinski curve. 












10 


HUI RAO AND SHU-QIN ZHANGf 


4. The chain condition, proof of Theorem 11.21 

Let (.4., r, 0, -<) be an ordered GIFS. Denote the invariant sets by For an edge 

a; € r, recall that is the associated similitude and t(a;) is the terminate state. 

For z G .4., a path u G F“ is called the lowest path, if is the lowest path in F” for 
all n; in this case, we call a = vrj(cj) the head of Ei. Similarly, we define the highest path 
cj' of Ff^, and we call b = vrj(ca') the the tail of E^. 


Definition 4.1. An ordered GIFS is said to satisfy the chain condition, if for any i € A, 
and any two adjacent edges ti ;,7 G F* with a; 7 , 

ff^(tail of = 07 ( head of £^*( 7 )). 


Glearly, for an ordered IFS {5i,..., S^}, the lowest coding is 1°° and the highest coding 
is N°°. Therefore, the head of K is the fixed point of 5i, denoted by Eix{Si), and the 
tail of K is Eix{S]sf). Consequently, the chain condition hold if and only if 


(4.1) 5i+i(Fix(5i)) = 5i(Fix(S'Ar)) for z = 1, 2,..., - 1, 

Condition (14.11) first appeared in Hata [15] , when dealing with the 1 /s-Holder continuous 
parametrization of self-similar sets. 


Proof of Theorem \l.Si Suppose {A,T,Q,A) satisfies the chain condition. Let z G A, 
and let u and 7 be two adjacent paths in F” with a; ^ 7. Let 77 = a; A 7 be the largest 
common prefix of cj and 7. Then u and 7 can be written as a; = rjojk+i • • • Wn and 
7 = ?77fc+i ... 7n- The fact cj and 7 are adjacent implies that 

(i) oJk+i and 7a;+i are adjacent edges in Fj where j = t{r}), and oJk+i ■< 7fc+i; 

(ii) ijJk +2 ■ ■ - OJn is the highest path in and 7^+2 • • • 7n is the lowest path in 

■pn—k—1 

^ hik+i)' 

By item (ii), we have h = (tail of G since the coding of b is initialled 

by iWk +2 Hence 

5 ajj,+i(tail of G 


Similarly, 


57,^1 (head of £4(7,+!)) ^ £7^+1...7„ ■ 


Therefore H £ 7 ^+ 1 ... 7 „ / 0 by the chain condition. So 

Eu> n £7 = (7T7(£tjfc+i...aj„ n £ 7 j._|_j^,, 7 ^) a 0) 


which proves that the GIFS is linear. 
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On the other hand, assume that (^, F, G, -<) is a linear GIFS. Fix i € A. Let loi and 71 
be adjacent edges in Fj satisfying cji -< 71 . Let {u}k )'^^2 highest path in F“^^^ and 

{lk )^=2 lowest path in Denote 

u\k= 0 Ji... LOk and 7 ]^ = 71 ... 7 ^, 

then for all k > 1, u\k and 7 !^ are adjacent path in F^ and so n / 0. As we 
know that 

(tail of = 7r^((wp)p>i) G 

5^1 (head of £'^(.^ 7 ) = 7 r„(( 7 p)p>i) G E^\^ 

for all A: > 1 , so the distance between (/^^^(tail of L't(a;i)) and ^..^^(head of £^*( 77 ) can be 
arbitrarily small. Thus 

(tail of £^t(a;i)) = 5-71 (head of £^*(^ 7 ), 

and the chain condition is verified. The theorem is proved. □ 

Corollary 4.2. An ordered IFS {S'!,..., S'tv} is a linear IFS if and only if (j4.1l) holds. 

5. Path-on-lattice IFS on the plane 

In this section, we study the path-on-lattice IFS on the plane (which we denote by C). 
Let L = Z -|- iZ be the square lattice or L = Z -|- wZ be the triangle lattice in the plane, 
where u) = exp(27ri/3). We define two points in L to be neighbors if their distance is 1. 
Then we obtain a graph and we still denote it by L. 

Let P be a path in L passing through the points 0 = zq; ^i) ■ ■ ■ j ^n-ij Zn = d \n turn. 
Let <I> = {fi>k}'k=i hs an ordered IFS on C such that 

(5.1) 0fc({O,d}) = { 2 ;fc_i,Zfc}, for all A: = l,...,n. 

We call such $ a path-on-lattice IFS with respect to the path P. Clearly the mapping (pk 
has the form (pkiz) = az + (3, or cpkiz) = az + ft with a, /? G C, and there are four choices 
of (pk for each k. If we indicate the four mappings by line segments with a half-arrow, 
then the IFS can be described by a path consisting of marked line segments. If all pk are 
of the form az -\- ft, then we say $ is reflection-free. 

Theorem 5.1. The path-on-lattice IFS is either a linear IFS, or its invariant set can he 
generated by a linear GIFS with two states. Moreover, the linear GIFS satisfies the OSC 
if the original IFS does. 
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Figure 7. Paths for Heighway dragon curve, Peano curve and Hilbert curve. 


Proof. Let {(l)k}k=i be a path-on-lattice IFS defined by (I5.ip . Let K be the invariant set. 
For k = 1,... ,n, define 


(5.2) 


Vk 


1, a {(j)k{0),f'k{d)) = {zk-i,Zk), 
-1, {4>k{0),4>kid)) = izk,Zk-i). 


We define an ordered GIFS with two state {1,-1} as follows: 


El — 4>l{Evi) + • • • + 4>n{Evn): 

E-l = (fniE-Vn) + • • • + (flliE-vi). 


(One may think that the first equation is corresponding the path P, and the second 
equation is corresponding to the reverse path of P.) Clearly Ei = E-i = K. 

We shall show that the head and tail of Ei are 0 and d respectively, and the head and 
tail of E-l are d and 0 respectively; then using (15.2p . we deduce that the GIFS (j5.3p 
satisfies the chain condition. According to ui = ±1 and = ±1, we have four choices. 

Let us denote the k-th edge emanating from Ei by Xi^k- If ui = 1 and = 1, then 
the first edge emanating from the vertex 1 is a self-edge, and hence (Ai^i)°° is the lowest 
coding. It follows that the head of Ei is 0. Similarly, the highest coding emanating from 
vertex 1 is (Ai^n)°°, and so that the tail of Ei is d. By the same argument, the head of 
E-l is d and the tail of E-i is 0. 

The other three cases can be proved in the same manner. Moreover, if all Vk equal 1, 
or all Vk equal —1, then GIFS (15.31) degenerates to a linear IFS. 

As for the open set condition, if the original IFS satisfies the OSC with an open set U, 
then GIFS (15.3p satisfies the OSC with open sets {U, U}. The theorem is proved. □ 


Remark 5.1. Algorithms of checking the open set condition of path-on-lattice IFS are 
discussed in [28]. For the path-on-lattice IFS (15.ip . if n = ||(i|p and the OSC holds, then 
K is a reptile (see Remark EH). 
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Example 5.2. We give several space-filling curves generated by path-on-lattice IFS’. All 
of them are reflection-free. The visualizations will be explained in next section. 

(1) Heighway dragon. The IFS is given by the path in Figure [7] (left), that is, 

^ + (1 + *)■ 

Denote the Heighway dragon by H. By Theorem 15.11 vi = 1 ,U 2 = —1, and {Ff, Ff} are 
the invariant sets of the following two-states linear GIFS: 

f El = (t>i{Ei) + ct)2{E-i) 

\ E-i = 4>2{Ei) + ())i(E_i), 

(2) Peano curve. The IFS is given by the path in Figure [T^middle) . Then vi = ■ ■ ■ = 

vg = 1, and it is a linear IFS: E = 4>i{E) + 4>2{E) (j)g{E). 

(3) Hilbert curve. The IFS is given by the path in Figure (T^right). Clearly ui = U 4 = 
— 1 ,U 2 = U 3 = 1, and the corresponding linear GIFS is: 

( El = 4>i{E_i) -I- 4>2{Ei) + 4>3{Ei) + (p4{E-i) 

\ E-i = (j)4{Ei) -|- (j)'i{E-i) -|- (/)2(E_i) -|- (j)i{Ei). 

Example 5.3. Gosper curve and anti-Gosper curve. The IFS of the Gosper curve is given 
by the path in Figure [ 8 ] (top-left). Clearly ui = U 4 = U 5 = ue = 1, V 2 = V 3 = vj = —1, and 
the corresponding linear GIFS can be obtained accordingly. 

If we forget the arrows of a path P, then we obtain a broken line and we call it the trace 
of P. It is shown in |28| that, among the path-on-lattice IFS’ which are reflection-free and 
have the same trace as the Gosper curve (there are 128 of them), none of them satisfies 
the open set condition except the Gosper curve and anti-Gosper curve. (The path of the 
anti-Gosper curve is determined by ui = U 2 = U 5 = 1 and V 3 = = vq = vj = —1. See 

F igure [HKbottom) .) 

6. Visualizations of space-filling curves 

We consider the visualizations of linear GIFS’ in this section. 

Let us start with linear IFS’. Let S = {Si}^i be a linear IFS satisfying the open 
set condition. According to Theorem 1.1, an optimal parametrlzation ip oi K can be 
constructed accordingly. To visualize the limit curve tp, we need to choose an initial 
pattern; indeed, the initial pattern can be any curve, but a suitable choice will make the 
visualization beautiful. 

Let us denote by Lg the initial pattern, and let a and b be its initial and terminate 


point, respectively. (Often Lq is chosen to be a line segment, which we denote by [a, 6].) 
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Figure 8 . Marked paths for Gosper island and anti-Gosper island. 

Fix an n > 1. For u S ,N}^, we denote by = Si^{x), and denote the 

follower of a; (if cj / (A^)")). Connecting Si^{a) and Si^{b) by Sui{Lq), and connecting 
Sui{b) and £'^^+( 0 ) by a line segment, we obtain the curve 

Ln — ^ ^ ( £a)(Fo) “1“ )• 

\u}\=n 

Here we use ^ to indicate that is the joining of small curves, where the order in given 
by uj. We call Ln the n-th approximation of the space-filling curve 99 . Different choices 
of initial patterns may give very different approximations in appearance, though the limit 
curve is the same. 

Example 6.1. Peano curve. The linear GIFS structure is given in Example 15.21 Figure 
[9Kleft and middle) shows two different 2nd approximations of the Peano curve. For the 
left picture, the initial pattern Lq = {3/2} is a singleton, and for the middle picture, 
Lq = [e, 3 — e] with e = 0.4. 

As for a linear GIFS, we have to choose an initial pattern for each Ej, which we denote 
hy Lq, , Lq . Denote the initial point of the pattern Lq by aj and the end point by bj. 
We define the n-th approximation of Ej to be 

^€r» 

where 6 ^ = and a^+ = gu,+ {at{Lj+))- 
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Figure 9. Left and middle: Visualizations of Peano curve. Right: Initial 
patterns of the Sierpihski curve. The red, green, purple and blue line 
segments are initial patterns for Ti,T 2 ,T 3 and T 4 , respectively. The dotted 
square is the unit square. 


Example 6.2. For space-hlling curves in this example, the linear GIFS structures are 
given in Example 13.5115.21 and 15.31 The visualizations and initial patterns are listed in the 
following tablq^. 


Hilbert curve 

Fig. UK middle) 

3rd 

-^0 — -^0 — + i} 

Heighway dragon 

Fig. ISKleft) 

00 

iJ = iS = {(l + i)/2} 

Gosper curve 

Fig. ISKright) 

3rd 

Lg = [ 0 , d], Lg = [d, 0 ] where d = 2 -|- 

Sierpihski curve 

Fig. (Upright) 

3rd 

Fig. fright) 


7. The four-star tile 

The four-tile star is a 4-reptile generated by the IFS 

Si{z) = - 2 ’ = “2 “ ^ “2 exp(i-^), S^{z) = “2 + exp(i-). 

Actually, the Sj map the big dotted triangle to the four small triangles in Figure 01 
m gave a visualization of the four-star tile without any explain (see Figure 0|). The 
mathematical theory behind are provided in [7], and here we give a sketch of it. 

7.1. Skeleton. We introduce the notion of s^e/eton of a self-similar set in [7] . (A skeleton 
is a kind of vertex set of a fractal.) It is shown that {ai,...,a6} is a skeleton of the 
four-star tile, where oi = |exp(i^), 02 = |i, 03 = 04 = 05 = ^6 = and 

uj = exp(i27r/3). (See Figure [10] and HU) 


^Figure [T] middle) shows the Hilbert curve generated by the path-on-lattice IFS multiplied by 1/2. 
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Figure 10. Another visualization of the four-star tile. 


7.2. Substitution rule. Let us consider the polygon P with vertices oi, 02 ,..., og. Here 
we regard P as a graph with directed edges, where each edge has the clock-wise orientation. 
Let us denote the edges by x, y, z, u, v, w. (See Figure fTTH eftl.l Figure [Upright) indicates 
the graph 


5i(P)u52(P)u53(P)u54(P) 

which consist of 24 directed edges. A Eulerian path of the graph is indicated by Figure 
[Upright), and the path is divided into 6 parts indicated by different colors. 

Replacing a segment in Figure [TlTlefti by the broken lines in Figure [TlT righti with the 
same color, we obtain a substitution rule. 



Figure 11. Substitution rule of the Four-star tile 
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7.3. Linear GIFS.. To precise the meaning of the substitution rule, we introduce the 
following linear GIFS: 

' X = 53(17) + 53 (F) + SsiW) + SiiU), 

Y = 5i(F) + Si{Z) + SiiU) + 54 (F), 

Z = 54(IF) + 54(X) + 54 (F) + 5i(fF), 

[/= 5i(X) + 52 (F) + 52(TF) + 52(X), 

F = 52 (F) + 52 (F) + S2iU) + 5i(F), 

^ W = SiiZ) + SsiX) + 53 (F) + 53 (F). 

It is shown that the four-star tile coincides with XUFUFUf/UFUlF, and it is a 
non-overlapping union (in Lebesgue measure). 

7.4. Visualizations. FieurefT^lefti provides the initial patterns of the third visualization 
in FigureHl If we choose the initial patterns in Figure fT^ rightl. we obtain the visualization 
in Figure [TOK right). 




Figure 12. Initial patterns for four-star tile. 


8. Proof of Theorem 11.11 Measure-recording GIFS 

In this section, we show that the invariant sets of a linear GIFS with the open set con¬ 
dition admit optimal parameterizations. An auxiliary GIFS, called the measure-recording 
GIFS, will play an important role. 

8.1. Preliminaries to dimensions and measures of graph-directed sets. Let (A, F, Q) 

be the GIFS given by (13.11) . 

We say a directed-graph (A, F) is strongly connected, if there exists a path from i to j 
for any pair i,j € A. 
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Let Te denote the contraction ratio of the similitude ge associated with e £ L. Define a 
matrix M(t), t > 0, as 

(8.1) M(t) = [ r* 

yeSFij 

Then there exists a unique positive number 6 , called the similarity dimension of the 
GIFS, such that 

p{M{5)) = 1 

where p{M) denotes the spectral radius of a matrix M. (See [20], [Ej.) 

Theorem 8.1. f|2nj ) Let (.4, T, Q) he a GIFS satisfying the OSC and let 6 be the similarity 
dimension. Then 

(i) dim// Ei = 5; and 0 < H^{Ei) < oo for all i = 1,2,... ,N ifV is strongly connected, 
(a) {H^{Ei), ... ,H^{E]\[)y is an eigenvector of M{5) corresponding to eigenvalue 1. 
(Hi) 'H^{Ei^ n E~^) = 0 for any incomparable cj ,7 G T*. (Two paths are said to be 
comparable if one of them is a prefix of the other.) 




In the rest of the section, we will always assume that Q satisfies the OSC, and that 
0 < H^{Ei) < oo for all i = 1,..., A^. Let us denote 

hi = TL^{Ei) and Pi = TL^\Ei, i = l,...,N. 

Now, we define Markov measures on the symbolic spaces T^, i ^ A. For an edge e G F 
such that e G Fjj, set 


( 8 . 2 ) 



s 

e • 


Using Theorem 18. Ilf hi. it is easy to verify that (pe)eGr satisfies 


(8.3) yy Pe = 1) for all i & A. 

We call (pe)eer a probability weight vector. Let Pj be a Borel measure on F“ satisfying 
the relations 


(8.4) Fi{[u}i . . . UVn]) = hiPoji .. . Puj„ 

for all cylinder [cui... ujn] ■ The existence of such measures are guaranteed by (18.311 . We 
call {Pi}(^i the Markov measures induced by the GIFS Q. The following result is folklore, 
see for instance [20l [18] . 
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Theorem 8.2. Suppose the GIFS Q satisfies the OSC and 0 < hi < +oo for all i. Let 
vTj ; —)• Ei he the projections defined by (13.3p . Then 

pLi=¥iO 7r“^ 

8.2. Measure-recording GIFS of a linear GIFS. Let {A,T,Q, be a linear GIFS 
such that the open set condition is fulfilled and 0 < hi = H^{Ei) < oo for all i. Set 

Fi = [0,hi], i = l,...,N. 

Fix a state i. We list the edges in Fj in the ascendent order with respect to 

71, ■ ■ ■ ,7£i- 

Recall that t{'y) denotes the terminate state of an edge 7. Then according to the set 
equation form of Ei can be written as 

Ei = H + 9■yt^ (Eti'Ye.))- 

Let 

fik ( 2 ^) = ^4^ + ^ ^ <£i 

be similitudes such that 

(8-5) Fi = /.yj -I -I- (-^ 4 ( 77 )) 

where the right hand side is a non-overlapping union of consecutive intervals from left to 
right. Indeed, we must have bk = i|8.5p holds by equation (18.3p . Doing 

this for all i ^ A, then p8.5p give us an ordered GIFS with the natural order. We denote 
this GIFS by 

and call it the measure-recording GIFS of (7l, F,^,-<). 

Glearly, the measure-recording GIFS inherits the graph structure and the order structure 
of the original GIFS; moreover, it records the Hausdorff measure information of the original 
GIFS. The following facts are obvious. 

• {Ei}^i are the invariant sets of the measure-recording GIFS. 

• For an edge e € F, the contraction ratio of fe is rf, and the similarity dimension 

of G* is 1. 

• G* satisfies the OSG. 

• The measure-recording GIFS shares the same symbolic spaces with the original 


GIFS. 
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Let 

TTi : rr ^ Si and Pi : r-^ Fi, i = l,...,N, 
be projections w.r.t. the GIFS (Q) and {G*), respectively. (See (13.3p .) Then 


Lemma 8.1. The Markov measure induced by the measure-recording GIFS coincides with 
that induced by the original GIFS. 


Proof. Let (pe)eGr and (p*)eer be the probability weights corresponding to G and G*, 
respectively. Since 


the two systems dehne the same probability weight vector and hence dehne the same 
Markov measure. □ 


Define 

(8.6) := TTj o pr^ 

The following lemma verihes that ijji is a well-dehned mapping from Sj to Sj. 



Lemma 8 . 2 . Suppose x E Fi has two pi-codings, say ^(x) = {0^,7}. Then 

7ri(cj) = 7ri(7). 

Proof. Write u = {uik)^^i and 7 = (7fc)^i. We claim that oji.. .ujn and 71 ... 7^ are 
adjacent for all n > 1, for otherwise, there exists pi ... such that 

Wi . . . Wn 7 Pi . . . Pn 7 71 . . . 7n, 

and the interval separates and contradicting to Pi(cj) = Pi(7) = x. 

Our claim is proved. It follows that n / 0, since (M, r,^,7) is a linear 

GIFS. Hence |'/rj(a;) — 7rj(7)| < diam + diam so the distance between 7ri(cj) and 
7rj(7) can be arbitrarily small, which implies that vrj(cj) = vrj( 7 ). □ 
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Now, we prove Theorem [LT] by showing that the mapping is an optimal parametriza- 
tion of Ei. 

8.3. Proof of Theorem II.11 Let Q* be the measure-recording GIFS of Q. Let ij^i = 
TTj o p~^. Let Vi = C\Fi be the restriction of the Lebesgue measure on Tj, pi = 'H^\Ei-, and 
Pj be the common Markov measure of Q and Q*. Then Pj = Pj o /ij = Pj o by 
Theorem 18.21 

(i) First, we prove that ifji is almost one to one. 

Let Qi be the set of points in Ei possessing more than one vTj-codings. Since 

Qi=[j U 

n>l 

and p,i{Ej n E^^) = 0 by Theorem 18. ll iiil. we obtain p,i{Qi) = 0. Denote 

then TTj is injective when restricted to F^ \ Aj, and Pi(Aj) = pi{Qi) = 0 . 

Similarly, let Q[ be the set of points in Ei possessing more than one pj-codings, then 
Vi{Q'j) = 0. Let A' = p~^{Q'i), then pi is injective when restricted to Ff^ \ A', and 
P,(A 9 = i.,(Q 0 = O. 

Let 

El = p,(rr \ (A, U A')), E' = vr,(rr \ (A, U A')), 

Then ijji : F' —>• S' is one-to-one, Vi{Fi \ F[) < Pi(Ai U A') = 0, and pi{Ei \ S') < 
P*(A,UA9 = 0. 

(ii) Secondly, we prove that "0^ is measure-preserving. For any Borel set B C we 
need to show that Vi{'ijj^^(B)) = Pi{B); due to (i), it suffices to show this hold for B C S'. 
Indeed, for B C E[, we have 

= Vi{pi o 7rr^(S)) = Pi o p-^ o Pi o 7rr^(S) = P* o 7r^"^(S) = pi(S), 

where the third equality holds since pi is a bijection when restricted to (F^ \ (Aj U A')). 
Similarly, for any Borel set S C M, one can show that /ii(i/’i(S)) = Vi[B). 

(iii) Finally, we prove the l/(5-Holder continuity of i/ii. 

Let xi,X 2 be two points in Si = [0,hi]. Let k be the smallest integer such that xi,X 2 
belong to two different cylinders of rank k, say, xi € pi([cj]), X 2 S Pi([7]), where cj / 7 G 
F^. It is seen that u = oji.. .ojk and 7 differ only at the last edge, that is. 


7 = Wl . ..bJk-llk- 
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We consider two cases according to u and 7 are adjacent or not. 
Case 1. u and 7 are not adjacent. (See figure [131) 


F F F 

^ (jj T] ^ •y 


Xi 


X2 


hi 


Figure 13. cj and 7 are not adjacent 

Then there is a cylinder rj = ui... Uk-iVk between u: and 7 , so 

— X 2 \> diam Fr/ > h ■ > h ■ rf,* • 


where u* = ui... Uk-i is the path obtained by deleting the last edge in lj, and 

h = min{hj; z = l,...,A^}, r^in = min{re; e G F}. 

Since xi,X 2 belong to / 9 j([cj*]), the images of xi and X 2 under vTj o which we denote 
by ui and y 2 respectively, belong to 7 rj([cj*]) = It follows that 

(8.7) \yi - y 2 \ < diam < D ■ r^* < D ■ ■ |xi - X 2 \^^^, 

where 

D = max diam Ei. 
l<i<N 

Case 2. u and 7 are adjacent. (See figure [HKleft).) 




Fv 


-I- 


Xi X3 X2 


Fui' Fi^n Fy 


Xi X3 X2 


Figure 14. u and 7 are adjacent 

Let X 3 be the intersection of and Ej. Let E be the smallest integer such that xi 
and X 3 belong to different cylinders of rank k', say, xi G pi{[u'\) and X 3 G pi{\u'']) (see 
Figure [CT rightll. then |xi — X 3 I > diam E^^ll since X 3 is an endpoint. 

Let 1/3 = Ei{x 3 ). Similar to Case 1, we have 

\yi - 2/3I < -D • ■ \xi - X3|^/^ 


By the same argument, we have 
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Hence, by the fact X 3 locates between xi and X 2 , 

(8.8) \yi - y 2 \ < 2D ■ ■ |a:i - 

Therefore, (|8.7p and (| 8 . 8 p verify the 1/6- Holder continuity ol i/i. □ 

Remark 8.3. We note that the initial point is the head of Ei, and the terminate 

point ipiihi) is the tail of Ei. 

Remark 8.4. In some text books, the space-filling curves is concerned; for example, ‘Real 
Analysis’ by Stein and Shakarchi (2005), ‘Topology’ by Munkres (2000), ‘Basic Topology’ 
by Armstrong (1997). The above proof extends the arguments in these books. 
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